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We investigate the effect induced by variations in the density profile of the Earth’s interior using 
a long-baseline neutrino oscillation experiment. 

At first, we point out two facts, (i) The most essential part of the matter profile is the hrst 
Fourier mode of the matter prohle function; and (ii) The Earth models based on the study of 
seismology include a large uncertainty for the first Fourier mode. Next, we show that there is a 
strong correlation between the average density value and the hrst Fourier coefficient in the analysis 
of oscillation probability. This means that the matter profile effect induces added uncertainty for 
the average matter parameter. 

Taking into account this extra uncertainty, we make numerical calculations for the sensitivity to 
CP violation search and show that CP sensitivity is impaired by this added uncertainty within a 
large U^s region. 
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I. INTRODUCTION 


The effects of the Earth’s matter Q play a very important role in long-baseline and high-energy oscillation experi¬ 
ments, such as a neutrino factory [||. Furthermore, the effect of variations in matter density and chemical composition 
along the baseline on the oscillation probability, which we call the matter profile effect, is also a controversial issue 
and has been debated in various contexts ||, ||, |^, |[ |[ ^ |^, Q |l^ . 

In the search for CP-violating effects, it is important to understand the structure of the Earth so as to estimate 
precisely the fake signal induced by the Earth’s matter. Some research |]^ has concluded that it is hard to obtain 
information about the interior of the Earth through neutrino experiments using the currently assumed size of statistics 
and realistic detector conhguration. Therefore, the science of geophysics has been applied to support the analysis. 
Much analysis has used the Preliminary Reference Earth Model (PREM) |l^, an Earth model which is based on 
seismology. When we deal with such a model, we need to be conscious of the fact that the model includes uncertainty. 
In this paper, we discuss the fact that this uncertainty reduces the sensitivity of the CP violation and we show that 
it is not a small effect. 

The procedure used to demonstrate our argument is as follows: Firstly, we recapitulate the method of the Fourier 
series expansion of the matter profile function 0 0 in Sec.^. Using this method, we show that only the first 
few modes are important in a high-energy experiment. Furthermore, we mention the uncertainty included in the 
seismological Earth models. We point out that a few percent error for the profile function can actually be interpreted 
as huge uncertainty for the Fourier coefficients, which can affect the CP sensitivity. Next, in Sec.HI, we show that 
there is a strong correlation between the constant matter parameter and the first Fourier coefficient of the matter 
prohle function within a wide energy and baseline region. This fact means that the uncertainty of the matter 
prohle effect brings added uncertainty to the average matter density. Then, we present numerical calculations to 
show quantitatively a loss of sensitivity induced by the uncertainty of the matter profile effect in Sec. 0 Finally, a 
summary is given in Sec0. 


II. METHOD OF THE FOURIER SERIES 

We would like to point out two facts in this section. (1) The most essential part of the matter profile effect is the 
first Fourier mode of the matter profile function. (2) The seismological Earth models include great uncertainty for the 
first Fourier mode. These facts force us to reconsider analyses in the baseline region which have so far assumed that 
the matter profile effect is not significant. In the case where the baseline length is 3,000 km, PREM tells us that the 
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matter profile function is almost flat and hence it can be expected that the matter profile effect is small. However, if 
we consider the uncertainty of PREM, how will the analysis change? 


A. Introduction of the method 


To see the effects induced by the matter profile, we derive the analytic expression using the method of the Fourier 
series UE- Expanding the matter profile function into the Fourier modes, we obtain an extremely clear viewpoint 
for the resonance conditions between the oscillation lengths of the neutrino and the matter profile undulation. By 
this expansion we can understand which modes, and what structures, are effective*. 

Now, we introduce our calculation method. We assume three generations and parameterize the mixing matrix of 
the lepton sector, Maki-Nakagawa-Sakata-Pontecorvo matrix (MNSP matrix), as follows; 


Uai — 
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{a = e,n,T, i = 1,2,3). (1) 


Here, s and c denote sin and cos, and HMajorana is the so-called Majorana phase matrix, which does not contribute to 
the oscillation phenomena. The effective Hamiltonian for neutrino propagation is 
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where is the squared-mass difference between the Hh and jth generations in a vacuum. We separate the matter 

effect into two parts, ao and Sa(x). The effect of the average matter is denoted by uq, and Sa(x) is the matter profile 
part, that is, the deviation part from the average density which depends on the position, x. After the separation, we 
expand the matter profile part into the Fourier series. 
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Note that the relation, a„ = al„ is always satisfied because of the condition that Sa{x) is real. For anti-neutrino, ag, 
a„ and i5 should be replaced by —Ug, —a„ and —(5, respectively. 

We treat Am^i and 6a{x) as perturbations and obtain the oscillation probabilities up to the first order of them, 
for example fl^. 
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The Fourier expanded-matter profile can not reproduce the boundary between two layers precisely. However, it will be shown below 
that the higher modes which construct the edge are not effective in high-energy experiments. 
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Here the effective mixing angle in the average density of the matter, </>, and the squared-mass eigenvalue, A±, in the 
average density are given by 
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Note that the effects of the asymmetric profile do not appear in the first order of the perturbations. 

The resonance conditions, A+ — A_ = 2Epn, show the energy range and width of the resonance induced by each 
Fourier mode. We understand that the higher Fourier mode can resonate only with the lower energy neutrinos whose 
oscillation lengths are shorter Q. Furthermore, the half-width of the amplitude becomes narrow as the mode becomes 
higher. This means that if the lower energy neutrino were to be observed precisely, the fine structure of the Earth 
could be known, although it is actually extremely difficult to achieve such observation. Therefore, we conclude that 
only the first few Fourier modes, which are determined by the large structure of the matter profile, are relevant in 
the currently assumed experimental set-ups. This is consistent with the results obtained using different methods 

O, iniTT 


B. Uncertainty of the Earth model 

Knowledge of geophysics is essential since it is very difficult to ascertain the profile of the Earth from a neutrino 
experiment. So far, PREM has been regarded as the absolute model. We have tended to expect that the error and 
the effect induced by the error are so small that they can be neglected without a careful consideration, although we 
have to estimate how much error the model includes in order to use the seismological Earth model. In order to discuss 
this error, we introduce another Earth model, akl35-f 0- 

Figure represents the matter profile function calculated using PREM (solid line) and akl35-f (dotted line) in 
the case where the baseline length is 3,000 km and 7,332 km, respectively. The profile function based on PREM is 
almost flat at L = 3,000 km, and this fact guarantees that the matter profile effect is small. According to akl35-f, 
however, it is not so flat. If we follow this model, we may be unable to ignore the matter profile effect, even at this 
baseline length. Furthermore, the authors of Ref. note that “the upper mantle density model should be treated 
with caution and may well change with further work.” The upper mantle and transition area (up to 670 km in depth) 
occupy a large part of the path of the neutrino beam in the case of L = 3,000 km. 

The baseline dependence of the average density and the hrst Fourier coefficient are compared in Fig.||. These two 
models are different in terms of the Fourier coefficient, although they are similar with regard to the average matter 
density. In particular, around L = 3,000 km, the difference of the first Fourier coefhcient is quite large. We need to 
recognize that seismological models include uncertainty, the size of which is at least equal to the difference in these 
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FIG. 1: The matter profile functions in the case of 3,000 km and 7,332 km which are calculated using PREM (solid line) and 
akl35-f (dotted line). 
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FIG. 2: The baseline dependence of the average matter density and the first Fourier coefficient which are calculated using 
PREM (solid line) and akl35-f (dotted line). Around L = 3,000 km, the average densities of these two models are almost the 
same within a few percent, but the first Fourier coefficients differ by more than 100%. 


two models^. The authors of Ref.[^ state that ±2 ~ 3% uncertainty in the upper mantle is a reasonable estimation, 
which produces a large, as much as 100%, uncertainty for the first Fourier coefficient. If 2.5% error for PREM is 
allowed, we can assume the matter profile depicted Fig.|| instead of PREM itself. This profile realizes 100% shift 
of the first Fourier coefficient without adjusting the average matter density. We would like to stress again that the 
uneertainty of the average density may indeed be small, but the uncertainty of the profile is not so small that it can 
be neglected, especially in L ^5,000 km. 


III. CORRELATION BETWEEN AVERAGE MATTER PARAMETER AND FIRST FOURIER 

COEFFICIENT 

Here, we will establish the existence of the correlation between the parameter for average matter density and the 
first Fourier coefficient of the matter prohle function. The effect induced by the matter profile and its uncertainty 
can be understood by the concept of this correlation. Although the average density can be determined at a few 
percent accuracy through studies of the seismic wave and gravitational effect etc., the matter profile effect adds extra 
uncertainty to the average matter density. As we pointed out in the previous section, the uncertainty of the first 
Fourier mode is not small, especially in the case where the neutrino beam passes mainly through the upper mantle 
and the transition zone, namely where the baseline length is around 3,000 km. 


A. One example: L = 7,332 km 

To corroborate that the correlation exists, we first present an example. In Fig.^, the oscillation probability for 
Ve —*■ is shown where the baseline length is 7,332 km. The solid line is calculated using the full-PREM profile. As 


t Assumption for the profile of the underground chemical component is made when the density profile is determined in geophysics. This 
assumption affects the density of the electron number that we really want to ascertain. 
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FIG. 3: An example where the deviation of the hrst Fourier coefficient from PREM is more than 100% in L = 3,000 km (dashed 
line). This is drawn by modifying PREM within 2.5%. The first Fourier coefficient for this profile is —0.084 g/cm®, which is 
almost twice as large as that of PREM, —0.043 g/cm®. Solid and dotted lines are the same as those in Fig.|^ 
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FIG. 4: Oscillation probability for I'e —> in L = 7, 332 km. The dotted, dashed, and dash-dot lines are calculated using the 
constant matter profile where the values of the matter parameter are 4.2, 4.4, and 4.6 g/cm®, respectively. The dash-dot-dot 
line is calculated by the up-to-first-mode profile where the values are po ~ 4.2 g/cm® and pi = —0.32 g/cm®. The solid line 
uses the full-PREM matter profile, and this line is quite similar to the dash-dot-dot and dashed lines. We set the oscillation 
parameters as sin o' = 0.5, sini/i = l/\/2, sin(fi = 0.1, Am§i = 2.5 x 10“®eV^, Am^i = 5.0 x 10“®eV^, and S = n/2. 


we showed in Ref. , this line is almost the same as the dash-dot-dot line which is calculated by the up-to-hrst-mode 
profile with po = 4.2 g/cm^ and pi = —0.32 g/cm^, whose values are based on PREM. 

The dotted, dashed, and dash-dot lines are those which are calculated using the constant-matter profile where the 
average density is 4.2, 4.4, and 4.6 g/cm^, respectively. The PREM tells us that the average density in this baseline 
length is 4.2 g/cm^. This figure shows that the constant prohle with the average density following PREM (dotted 
line) cannot reproduce the behavior of the oscillation probability with the full-PREM profile (solid line). However, 
Fig.Q also shows that if a shift in the constant value is allowed, then a good ht with the full-PREM calculation can 
be realized at po = 4.4 g/cvo?. This fact means that the shift of the constant matter parameter can copy the matter 
profile effect, which is almost equal to the effect induced by the first Fourier modes with this baseline length. We can 
expect that there is a strong correlation between the parameter for the average matter density and the first Fourier 
coefficient of the matter profile function. Indeed, this correlation is not an accidental phenomenon in this example. 

We now consider the mechanism on which this correlation is based. The existence of the correlation suggests that 
the common shift in the average matter parameter over a wide energy region can imitate the effect of the first Fourier 
mode. This statement means that the relation. 
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is satisfied by the constant Auq over a wide energy region. In the above example, the constant shift Apo = 0.2 g/cvo? 
works well within the energy region above 4 GeV. Since we can assume that the unperturbed term of the oscillation 
probability is dominant over the other perturbative terms and that the matter profile effect can be represented by the 
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It is not inconsequential whether or not the constant shift Aoq can maintain the relation. To clarify this issue, we 
divide the energy region into two regions where L ~ 0 and ^ IL are satisfied, and we investigate each 

region. 



FIG. 5: Energy dependence of for neutrino (solid line) and anti-neutrino (dashed line) in L = 7,332 km. The oscillation 

parameters used in this figure are the same as those used in Figg 

Figure ^ shows the energy dependence of corresponding to the case of Figg. In the region of ^+ 4 / L ~ 0 , 

the condition eq. ( 0 ) is approximated to 

Aao = -^Re[ai] + O | . (11) 

Using the value p\ = —0.32 g/cm^, we understand from eq. ( 0 ) that shifting po by 0 . 2 g/cm^ makes the oscillation 
probability calculated with the constant profile clearly mimic the real oscillation probability, including the matter 
profile. Within the A^^^A^p ^ E region, the condition eq.(p^ is approximated by 

Aao = -^Re[ai] + O . (12) 

This means that the shift of gq so as to mimic the matter profile effect is also about 0.2 g/cm^, which is the same as the 
shift in the L ~ 0 region. Therefore, in a wide energy region, such as both the AlzA^p ^ q and AtzA^p ^ E 

regions, the common shift of oq can indeed copy the effect of the first Fourier mode^. 

It is an essential point for the existence of the correlation that an experiment is only sensitive to these two energy 
ranges^. We note that the condition for anti-neutrino is the same as that for neutrino. Therefore, even if the analysis 
is made using both neutrino and anti-neutrino, the correlation still exists if the observed energy region for both 
neutrino and anti-neutrino satisfies either Atz^p ~ 0 or Atz^p ~ In the case of L = 7,332 km, the energy 
range for anti-neutrino is not around either AtzA^p ~ 0 or AtzA^p ~ so if we could observe an adequate number 
of anti-neutrino events, the correlation would cease to hold, although the anti-neutrino event can be expected to be 
too short to be significant in a statistical sense. Therefore, we can conclude that the correlation still exists for this 
baseline length. 


t Moreover, we can see that this common shift works well in the intermediate region by expanding the oscillation probability around 



§ In the region where —L ~ tt holds, the relation becomes Aag = — Re[ai] -|- O (^At^A—L — Aj , and hence the required shift to 

mimic the matter profile effect is significantly different. 
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B. Baseline region where the correlation exists 


There is a strong correlation between oq and ai within the energy and baseline region where 
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satisfied, as we have established in the previous sub-section. The energy and baseline dependence of 
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FIG. 6: Energy and baseline dependence of L for neutrino (left plot) and anti-neutrino (right plot). The reference values 

of the parameters are the same as in Fig.| 

neutrino and anti-neutrino is shown in Fig.o. This figure tells us the region where the correlation exists. In T < 5,000 
km, —is near 0 or ^ for both neutrino and anti-neutrino throughout almost all the energy range. Therefore, 
the common shift of Op for neutrino and anti-neutrino can mimic the effect of Oi, that is, the correlation does exist. 
In 5,000 km ^ L < 7,500 km, L for anti-neutrino is not around either 0 or However, since there will not 

be a significant number of anti-neutrino events, in this region, the statistics will be dominated by neutrino events 
and hence the correlation still exists. In the region beyond 7,500 km, the high-energy neutrino no longer follows the 
condition, and so the correlation will cease to hold. We conclude that there is a strong correlation between op and ai 
in the case where the baseline length is less than 7,500 km. 

We would like to note one more thing. According to eqs.(H) and (^), the same relation also holds for i/g Vr and 
Ve So, even if we utilize these channels, the correlation will not be broken. 


IV. NUMERICAL CALCULATIONS 


To illustrate how the uncertainty of ai impairs sensitivity to the CP violation, we present the numerical results for 
it, including the matter profile effect and its uncertainty. In an effort to clarify the difference from former research, 
we show simultaneously the result without consideration of the matter profile. 

First, we explain briefly the procedure for drawing the sensitivity plot: We define the test statistics. 


= min 
osc. param. 

(5 = {0,7r} 



jjyth X X NfY 

(7V‘h)2 X + (Vth)2 X 


(13) 


Here, iV®’' is the “expected number of events” for i/f, calculated using the full-PREM matter profile with 5 = 7r/2, 
iV*** is the “theoretical number of events”, calculated with the constant matter profile or the up-to-first-mode profile 
with 5 = {0,7r}, and N denotes the number of events for Df. —> The index i stands for the energy bin. The 

parameters contained in and iV*** are varied within given ambiguities. We adjust them and minimize to 
introduce the effect of the parameter correlation §0- The widths of uncertainty of the parameters concerning the 
atmospheric and solar neutrino experiments are expected to be narrowed by near-future experiments [po[ . Therefore, 
we assume 


A(sin^/;) = 1%, A(Am 3 i) = 3%, 

A(sina;) = 5%, A(Am 2 i) = 5%, (14) 

and for the other parameters including the matter effect we assume some options and compare them to each other. 
We do not deal with systematic error. We require 

> xi 9 %(d.o.f = bin) 


(15) 

















in order to claim that the hypothesis with i5 = {0, tt} is excluded at the 99% level of significance, where the right-hand 
side is the distribution function whose degree of freedom is the number of energy bins, not the number of parameters 
which fit. This is because we adopt the concept, the power of test (see appendix of Ref. fl^ for more details). The 
reason why we use this, less familiar, concept is that we firmly believe that we must pay attention to the fact that the 
best-fit point suggested by an experiment is not always located on the point chosen by nature. Remind yourself how 
the best-fit point for the solar neutrino deficit has changed. We actually know that there were several good-fitting 
regions for the solar neutrino experiments which were separated from each other on the parameter plane and which 
were not distributed only around the best-fit point. Moreover, the best-fit point itself moved from one region to 
another. To discuss the feasibility of observing some quantity with an experiment, we have to consider this fact. 

We also note in passing that the degree of freedom for so-called Ay^, which is often used in an estimation of the 
oscillation parameters, is not the number of the parameters since, for example, two strongly correlated parameters are 
not independent of each other, and hence we should count them as one parameter. Indeed, as is commonly known, 
there are strong correlations for some parameters in the case we deal with^. We should be more careful about what is 
a truly measurable quantity pl| , and which quantities statistics are sensitive to fl^ . For instance, if the statistics can 
be constructed so as to be sensitive to the CP violating effect, then the parameter is only one, Jarlskog’s parameter, 
and the degree of freedom should be one. Strictly speaking, only when the statistics are linearly dependent on the 
parameters, will the degree of freedom coincide with the number of the parameters. 


Amli[eV‘^] 



FIG. 7: Sensitivity reach for CP violation with 5 = n/2, a baseline length of 3,000 km, and a muon energy of 30 GeV. The 
theoretical number of events, of the dotted curve is calculated using the up-to-first-mode profile, where the widths of 
parameter uncertainty are assumed to be A(ao) = 3%, A(ai) = 200%, and sim^ = 0 ~ 0.16. For the other parameters, the 
uncertainty in eq.([l^ is assumed. The solid, dash-dot and dashed curves are calculated using the constant matter profile 
whose uncertainties are Aao = 0%,3% and 5%, respectively, and the uncertainty for the other parameters is the same as for 
the dotted line. 


To draw a sensitivity plot, we assume 40 kt detector and 10^^ muon decays in the neutrino factory scheme, and 
interpret eq.(p^ to the condition for Airi^i and sini)). Figure ^ indicates the lower boundary of Aml^ and |C43| 
(= sincj)) to reject the hypotheses that (5 is 0 and tt at a 99% level of significance, when nature adopts the value 
(5 = 7r/2 in the case where the baseline length is 3,000 km, muon energy is 30 GeV, the detection threshold is 5 GeV 
and the width of energy bin is 2.5 GeV (10 bin). The reference values are the same as those in Fig.^, except for Am^i 
and sinf] 

sin= 1/v^, Amg^ = 2.5 X 10~^eV^, sinw = 0.5, i5 = 7r/2. (16) 

The dotted curve is calculated whilst considering the matter profile effect and its uncertainty up to the first Fourier 
mode in the calculation of We take 3% for ag and 200% for oi as the widths of uncertainty. The other curves are 
calculated assuming the constant matter profile with a different uncertainty for ag. The solid, dash-dot, and dashed 
curves correspond to A(ao) = 0%, 3% and 5%, respectively. In the minimization process, sin^ in V*** is taken as an 
arbitrary value between 0 ~ 0.16, and the uncertainty of the other parameters is assumed, as in eq. O)- 

The dotted curve differs from the dash-dot curve in which the same uncertainty for ag is assumed, but that for oi 
is not included. In contrast, it is very similar to the dashed curve whose uncertainty for ag is larger (5%), but where 
a constant matter profile is assumed. This fact can be explained by the correlation between ag and oi eqs. (© and 
(|l^): Uncertainty of ai is translated into that of ag through the correlation, and this extra uncertainty gives rise to 
an extra absorption of the signal of the CP violating effect. We know that the effect of the fake GP signal induced by 


^ Depending on the conditions, it breaks the correlation to regard the oscillation parameters as not free ones, but as restricted ones. In 
the numerical calculation, this effect is automatically introduced by setting the widths of the uncertainty. 
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the matter effect becomes more serious as |C/e 3 | increases, so if \Ue^\ is measured just below the current boundary, we 
will have to take the effect induced by the density profile of the Earth more seriously. 

In the case of L = 3,000 km, the determination of sini^ and the reduction of its uncertainty do not contribute 
towards improving the sensitivity much, since the matter effect itself is large at this baseline length. Furthermore, 
lowering the detection threshold also hardly helps to improve the sensitivity because the number of events in the 
high-energy region is overwhelmingly greater than that in the low-energy region. If we can realize a lower detection 
threshold, an advantage can be gained with a shorter baseline length and a lower energy beam. 


Am2i[eV^] 
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FIG. 8: Sensitivity reach for CP violation when 5 = tv 12 in the case where the baseline length is 1,000 km, the muon energy is 
11 GeV, and the detection threshold is 1 GeV. The up-to-£rst-mode profile is adopted for the dotted curve where its widths of 
uncertainty are A(ao) = 3%, A(ai) = 200%, A(sin (f)) = 10% and eq.(|l4|) for the other parameters. The solid curve is calculated 
using the constant matter profile with A(ao) = 3%. By comparison with Fig0. it is obvious that the uncertainty of the matter 
effect is not serious in this situation. The dash-dot curve is calculated in a similar manner to the dotted curve, except that 
sincj) varies from 0 to 0.16. The dotted and dash-dot curves show that the uncertainty of sin0 will play an important role in 
this situation if a large sin 0 is established. 


Figure || shows the sensitivity in the case where the baseline length is 1,000 km, the muon energy is 11 GeV and the 
energy threshold is 1 GeV. In this plot, the dotted curve represents the calculation using the up-to-first-mode profile 
and its widths of parameter uncertainty are A(ao) = 3%, A(ai) = 200%, A(sin0) = 10%, and eq. ©) for the others, 
the solid curve is calculated by assuming a constant matter profile with A(ao) = 3%, and the dash-dot curve is the 
same as the dotted one, except that sincj) in is assumed to be an arbitrary value between 0 ~ 0.16. 

Within the large \Ue 3 \ region, the sensitivity is better than that in Fig.^ Furthermore, the sensitivity of this choice 
is more robust against the uncertainty of the matter effect than in the case of L = 3,000 km. In contrast, although the 
signal of the GP violation is also small in the small |C43| region, the fake GP violation effect induced by the matter 
effect is suppressed more strongly than the signal of the genuine CP violation, and hence, adopting a longer baseline 
and higher energy option may be advantageous to the CP violation search. 

This behavior does not come from the property of the statistics as defined by eq. When the uncertainty of the 
constant matter parameter is assumed to be larger, the optimization, using so-called Ay^, also suggests that a shorter 
baseline length and lower energy is better (see, for example O. Yasuda, in Ref|^). The existence of this correlation 
tells us that even though it is said that the uncertainty of the average matter density on the baseline is well estimated, 
this is not the entire uncertainty of the constant matter parameter. We would like to stress that the uncertainty of 
the matter effect is no longer so small when we introduce the matter profile effect. Therefore, we should regard it 
more seriously whatever statistics we use, especially in the case where the main part of the neutrino beam path is 
the upper mantle and transition zone, including a large uncertainty for the density profile. We need to deal with the 
matter effect much more cautiously. 


V. SUMMARY AND CONCLUSION 

We pointed out that there was a very strong correlation between the constant matter parameter, oq, and the first 
Fourier coefficient of the matter profile function, oi, within a wide energy and baseline region. This fact means that 
the uncertainty of ai is translated into that of ag, and this gives added uncertainty to ag. 

We also showed that there is a huge uncertainty in ai. This is due to the fact that seismological Earth models 
include a large uncertainty for the density profile in the region, the upper mantle and the transition zone, which is 
24 ~ 670 km in depth, and this uncertainty can cause a huge, even a few hundred percent, uncertainty for ai. The 
existence of the correlation suggests that this huge uncertainty affects CP sensitivity since, due to the correlation, the 
uncertainty in oi gives added uncertainty to og. 
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We present the sensitivity plot for the CP violation effect including the matter profile effect and its uncertainty. 
In the case where the baseline length is 3,000 km, most of the path of the neutrino beam is occupied by the upper 
mantle, which includes large uncertainty. We show numerically that 200% uncertainty for ai can be interpreted as 
about 2% extra uncertainty for oq, which should be added to the original uncertainty of oq, and this result confirms 
our expectations from the correlation. This extra uncertainty makes the CP sensitivity worse, especially within the 
large | C/e 3 1 region. A shorter baseline and lower energy option can avoid this disadvantage if the detection threshold 
can be lowered. On the contrary, if small |C/e 3 | is established, then a long-baseline and high-energy option may be 
better than a shorter baseline and lower energy, because of the statistics. 

We made some comments in answer to questions about the statistics which we used. We consider the fact that 
the best-fit parameter suggested by the experiments is not always distributed only around the parameter chosen by 
nature. Therefore, to discuss the feasibility of observing the CP violation effect, we need to consider this fact. This 
led us to use the concept, the power of test. 

In this study, we do not consider systematic uncertainty. Of course, in order to optimize the experimental config¬ 
urations, it is necessary to take account of systematic error. However, in any case, we can conclude that we should 
regard the uncertainty of the Earth’s matter as a more severe problem than that has so far been assumed. We need 
to be much more conservative when estimating the matter effect. 
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